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Abstract. We calculate the temperature and polarization patterns generated 
in anisotropic cosmological models drawn from the Bianchi classification. We 
show that localized features in the temperature pattern, perhaps similar to the 
cold spot observed in the Wilkinson Microwave Anisotropy Probe (WMAP) data, 
can be generated in models with negative spatial curvature, i.e. Bianchi types 
V and VII|j. Both these models also generate coherent polarization patterns. 
In Bianchi VII^, however, rotation of the polarization angle as light propagates 
along geodesies can convert E modes into B modes but in Bianchi V this is not 
necessarily the case. It is therefore possible, at least in principle, to generate 
localized temperature features without violating existing observational constraints 
on the odd-parity component of the cosmic microwave background polarization. 
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1. Introduction 

Observations of the temperature anisotropies of the cosmic microwave background, 
particularly those from the Wilkinson Microwave Anisotropy Probe (WMAP) [HE], 
form the foundations of the current ("concordance") cosmological model [3]. However, 
WMAP has also uncovered tantalizing evidence of departures from the standard 
framework. In particular, detailed analysis of the pattern of temperature fluctuations 
has led to the identification of a Cold Spot gHH H H H [9l [10] . The level of departure 
from isotropy is relatively small but it is highly significant from a statistical point of 
view. The presence of this feature seems to be inconsistent with the assumption of 
global isotropy upon which the concordance cosmology is based and could be evidence 
that we live in a globally anisotropic Universe, i.e. one not described by a Friedmann- 
Robertson- Walker (FRW) model. 

The Bianchi classification arranges all possible spatially homogeneous but 
anisotropic relativistic cosmological models into types depending on the symmetry 
properties of their spatial hypersurfaces [HI [12] . It has been known for some time 
that localized features in the radiation background can occur in Bianchi models with 
negative spatial curvature [131 [HI [13 [lEl [H] ■ The physical origin of such features lies 
in the focussing effect of spatial curvature on the geodesies that squeezes the pattern of 
the anisotropic radiation field into a small region of the sky. Only a few of the Bianchi 
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types contain the FRW model as a limiting case and, from this subset the model which 
appears to best able to reproduce the anomalous cold spot is the Bianchi VII/i case 
[ini HOI 1211 123] • However, as well as forming distinctive features in the temperature 
pattern, anisotropic cosmological models also generate characteristic signatures in 
the polarized component of the background radiation. Thomson scattering generates 
polarization as long as there is a quadrupole anisotropy in the temperature field of 
the radiation incident upon the scattering particle. In the concordance cosmology 
the temperature and polarization patterns are (correlated) stochastic fields arising 
from their common source in scalar and tensor perturbations arising from inflation. 
In a Bianchi cosmology, however, the patterns are coherent and have a deterministic 
relationship to one another owing to their common geometric origin. It has recently 
been shown [24l[25l|26] that the properties of the polarization field produced in Bianchi 
Vllft are inconsistent with the latest available WMAP polarization data [27] because 
they inevitably involve a large odd-parity (B-mode) contribution that exceeds the 
experimental upper limit. 

In a forthcoming paper we present an exhaustive study of the temperature and 
polarization anisotropics produced by those Bianchi types that possess an FRW limit 
[26] . The purpose of studying these models is to characterize as fully as possible 
the radiation fields they can produce in order to separate them as clearly as possible 
from residual systematics and thus provide the strongest possible constraints on exotic 
cosmologies. In the present paper, however, we focus on the very specific question of 
whether localized spots necessarily involve a large B-mode polarization. 



2. Bianchi Cosmologies 

The models we consider are based on Einstein's general theory of relativity and we 
use the field equations in the form 

Gab = Rab — -^Rgab = ^ab " ^9ab, (1) 

with Rab being the Ricci Tensor, R the Ricci scalar. Tab the energy-momentum tensor 
and A the cosmological constant. Indices a and b run from to 3. We use units where 
SttG — c — 1. In terms of a coordinate system x", the metric gab is written 

ds^ — gabdx°^dx^ = {hab — UaUb)dx°'dx^ , (2) 

where is the fluid velocity; the signature of gab is ( — h ++)• Starting from a local 
coordinate system x^ we construct a tetrad basis [H] 

e. = eL^ (3) 

such that 

9ab = e\elg,j = e\eu = diag(-l, +1, +1, +1) (4) 
meaning that the tetrad basis Ba is orthonormal. The Ricci rotation coefficients, 

Tahc = ^a^ci-J^i, (5) 

are the tetrad components of the Christoffel symbols; semicolons denote covariant 
derivatives. In general, the operators defined by equation ([3]) do not commute: they 
generate a set of relations of the form 

[ea^ebj^-fab^c. (6) 
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The Ricci rotation coefficients are just 

Tafcc ~ 2^1abc + Icab — "fbca) ■ (7) 

The matter flow is described in terms of the expansion a-nd shear Gab'- 

Ua-b = Wafc + -dab - UaUb 

(Jab = t^ab - ^-ffhab = l^ab ^ Hhab, (8) 

where d = Tr{'dab) = ^aa and the magnitude of the shear is = a^^Uab /2. We now 
take the time-hke vector in our basis to be the fluid flow velocity so that u°- — 6^ and 
Ua = — (JJJ. The remaining space-hke vectors form an orthonormal triad, with a set 
of commutation relations like that shown in equation ([6]) but with an explicit time 
dependence in the "structure constants" describing the spatial sections: 

[e.,e,] =7,^-(t)efc. (9) 

Without loss of generality we can write 

-f'[^=e,,in"' + S^a,-S'[a,, (fO) 

for some tensor Uij and some vector a^. The Jacobi identities require that UijU^ = so 
we choose — (a, 0, 0) and Uij — diag(ni, n2, "-s)- The four remaining free parameters 
are used to construct the Bianchi classification described in detail elsewhere [TTl [12] . 

We are interested in cosmological models that are close to the homogeneous and 
isotropic FRW case, but not all the Bianchi types contain this solution. Those that 
do are types I, V, VIIq, Vllh and IX. Bianchi I and Bianchi VIIq are spatially flat, 
Bianchi IX is positively curved and Bianchi types V and VII/j have negative spatial 
curvature. The open cases are of particular interest in this paper as they permit the 
focussing of anisotropic patterns into small regions of the sky. The scalar curvature 
R of the spatial sections is given in terms of the Bianchi parameters as 

R = — (2ni?i2 + 2rtin3 + 2^2713 — n\ — n\ — n\) — 6a^. (11) 

For Bianchi V we have ni = n2 = 713 = so that R = —6a^. In Bianchi Vllh we have 
ni = but 712 7^ and 77.3 7^ 0; the parameter h is defined hy h = 0^/772773. If we take 
772 = 773 = 77, the canonical form of the model, we also have R = — 6a^ = —6n^h. 

The models we consider have a single preferred axis of symmetry. The alignment 
of the shear eigenvectors relative to this preferred axis determines not only the 
dynamical evolution of the model through the field equations, but also the temperature 
and polarization pattern that gets imprinted into the cosmic background radiation. 

3. Temperature and Polarization Anisotropies 

The transfer equation for polarized radiation propagating through space-time can be 
described by a combination of propagation vector and a polarization vector. In the 
tetrad frame these can be written as 

in terms of the standard Stokes parameters, /, Q, U and V where the superscripts 
indicate the parts of the (complex) photon distribution involving spin weights 
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-— = e^-jk-n^^kik^ + a^hk^ - a^, (14) 



and 2. The direction of a light ray in three-dimensional space is defined by fc* = 
(cos 0, sin cos 0, sin sin 0) and we have also used the complex unit vector 

1 \ 8 id 

m^^-—dk^^- — {-^ + — -^)k\ (13) 

Neglecting shear-dependent terms in the geodesic equation, the following equation can 
be obtained [TH [15] that describes the evolution of the direction of propagation: 

1 dki 

£ dt 

£ is the photon energy. This leads to 
1 d9 

-— = [a + (^3 — 712) cos0sin6'] sin6' 

= [rii — + (ns — ^2) cos^ 0] cos 9. (15) 

£ dt 

The radiation distribution described by is affected by the properties of the 
geodesies along which photons propagate but is also altered by scattering. If only 
elastic scattering is included, this is governed by a Boltzmann equation 

VN^ = T(-iV^ + J^), (16) 

involving the operator 

I? = — - £7"— + ^(m*8 + m'S) + iS^T'^i^k^k'^e'''"', (17) 
ot 0£ -y/2 

where r the optical depth; the amplitude of the polarization depends strongly on t 
but we shall not explore its dependence in detail in this paper. The horizontal bar 
denotes complex conjugation. We also have j°- — T^-gfc* -I- r°--jk^y , so the non-zero 
connection terms are given by 

— — ?9 ■ 

«i ~ jo ~ "I'J 

r'jfc = -r^,-, = \iemn\ - euku] + e^jmi) + {S'^a^ - S}ak) (18) 

if we assume no vorticity (w'^ = 0) and no acceleration (F^Qj — iii). The Boltzmann 
operator pT|) involves terms that describe the propagation of the radiation field along 
the geodesies, the effect of the shear on the photon energy distribution, and the 
rotation of the polarization angle x defined by 2x — arctan(t//Q). The emission 
term J a (which involves terms describing the fluid flow) contains only harmonics up 
to Z = 2, since the radiation modes with I < 2 are damped as well as re-radiated by 
Thomson scattering, while higher-order modes / > 2 are only damped. 

It has become conventional to decompose the polarized component of the cosmic 
microwave background into modes classified according to their parity. The even modes 
are called E-modes and the odd modes are the B-modes. The latter are of particular 
interest in the context of inflationary cosmology as they cannot be sourced by scalar 
perturbations and are therefore generally supposed to be a signature of the presence 
of primordial tensor perturbations, i.e. gravitational waves [28l [29] 



= Nf^m'^ = (g - iU){n) = ^ a_2J™ -2i^™(ri), 



Ir. 



= Nf^m^i = (Q + iU){h) = 02 2Yimin) 



Im 
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Figure 1. Temperature and Polarization maps for Bianchi models with localized 
features (aligned with the known CMB Cold Spot). The left panels show the 
temperature, the right the total polarization P = ^/ _|_ (72^ The top row shows 
an example of Bianchi VIIj, with a compact spiral feature produced by focussing 
and twisting; the polarized component contains a significant B-mode. The bottom 
row shows a Bianchi V model in which the initial quadrupole is focussed but not 
twisted. The polarized field in the latter case has no B-mode. 



aE,/m — — (a2,im + a-2,im)/2, as,im = i(a2,im " a-2,im)/2 (19) 

The overall level of polarization (in both E and B modes) increases strongly with the 
optical depth to Thomson scattering r through 

However, the equations of the previous section for the evolution of the polarized 
radiation distribution allow us to establish some firm implications for the relative size 
E and B modes just by considering their symmetry. First, the scattering term J'^ 
produces a pure E-mode quadrupole anisotropy. However, depending on the initial 
conditions, the redshifting effect of shear can produce either E or B modes (just as 
a gravitational wave perturbation of FRW can). The physical processes described by 
Equations (14)-(17) can, in principle, convert E modes into B modes and vice versa. 
In Bianchi VH^ these effects are unavoidable so, even if there is no initial B-mode, one 
is inevitably generated as the Universe evolves. In Bianchi V, however, the last term 
in Equation (17) does not contribute to the mixing of E and B modes, at least at this 
level of perturbation theory. Since initial conditions exist in which the polarization is 
purely E mode there are therefore models of this type that can produce spots without 
any B modes. Figure 1, which was plotted using the Healpix software [301 1 shows 
illustrative examples of Bianchi Vllh and Bianchi V that show localized patterns with 
and without B-modes respectively. 

4. Discussion and Conclusions 

We have shown that it is possible, in a Bianchi V cosmological model, to generate a 
localized temperature anomaly qualitatively similar to the known Cold Spot without 
necessarily producing a large B-mode polarization. Observational limits on the B 
mode alone are therefore not sufficient to exclude global anisotropy as a possible 
explanation of the famous Cold Spot. Although the set of Bianchi models that can 
evade the limits set by B mode polarization is small, a more rigorous analysis requires 
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a fuller parametrization in which the Bianchi parameters are constrained with other 
cosmological parameters (such as the optical depth) using a complete description of 
the polarized radiation field. 

We stress that there are severe difficulties with anisotropic universes as 
explanations for the overall pattern of observed CMB anomalies. Most important 
among these is that a significant (negative) spatial curvature seems to be at odds with 
measurements that clearly prefer a flat universe [HE]. Nevertheless, these models may 
provide important clues that can lead to more effective and efficient use of observations 
to test exotic cosmologies. For example, there are other ways in which the polarization 
angle for CMB photons can be rotated [31]; the lack of any observed B-mode allows 
present observations to place strong constraints on such models also |32) . 
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